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Abstract 



There is a widespread recent interest in using ideas from statistical physics to model 
certain types of problems in economics and finance. The main idea is to derive the 
macroscopic behavior of the market from the random local interactions between agents. 
Our purpose is to present a general framework that encompasses a broad range of mod- 
els, by proving a law of large numbers and a central limit theorem for certain interacting 
particle systems with very general state spaces. To do this we draw inspiration from 
some work done in mathematical ecology and mathematical physics. The first result is 
proved for the system seen as a measure-valued process, while to prove the second one 
we will need to introduce a chain of embeddings of some abstract Banach and Hilbert 
spaces of test functions and prove that the fluctuations converge to the solution of a 
certain generalized Gaussian stochastic differential equation taking values in the dual 
of one of these spaces. 

1 Introduction 

We consider interacting particle systems of the following form. There is a fixed number 
N of particles, each one having a type w G W. The particles change their types via two 
mechanisms. The first one corresponds simply to transitions from one type to another at 
some given rate. The second one involves a direct interaction between particles: pairs of 
particles interact at a certain rate and acquire new types according to some given (random) 
rule. We will allow these rates to depend directly on the types of the particles involved and 
on the distribution of the whole population on the type space. 

Our purpose is to prove limit theorems, as the number of particles N goes to infinity, for 
the empirical random measures i/^ associated to these systems. is defined as follows: if 
77^ (j) e W denotes the type of the i-th particle at time t, then 
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where is the probabihty measure on W assigning mass 1 to w. 

Our first result, Theorem 1, provides a law of large numbers for on a finite time 
interval [0,T]: the empirical measures converge in distribution to a deterministic continuous 
path vt in the space of probability measures on W , whose evolution is described by a certain 
system of integro-differcntial equations. Theorem 2 analyzes the fluctuations of the finite 
system ly^ around i^t, and provides an appropriate central limit result: the fluctuations are 
of order 1/^/N, and the asymptotic behavior of the process y/N^v^^ — vt) has a Gaussian 
nature. This second result is, as could be expected, much more delicate than the first one. 

In recent years there has been an increasing interest in the use of interacting particle 
systems to model phenomena outside their original application to statistical physics, with 
special attention given to models in ecology, economics, and finance. Our model is specially 
suited for the last two types of problems, in particular because we have assumed a constant 
number of particles, which may represent agents in the economy or financial market (eco- 
logical problems, on the other hand, usually require including birth and death of particles). 
Particle systems were first used in this context in FoUmer (1974), and they have been used 
recently by many authors to analyze a variety of problems in economics and finance. The 
techniques that have been used are diverse, including, for instance, ideas taken from the Ising 
model in FoUmer (1974), the voter model in Giesecke and Weber (2004), the contact process 
in Huck and Kosfeld (2007), the theory of large deviations in Dai Pra, Runggaldier, Sartori, 
and Tolotti (2007), and the theory of queuing networks in Davis and Esparragoza- Rodriguez 
(2007) and Bayraktar, Horst, and Sircar (2007). 

Our original motivation for this work comes precisely from financial modeling. It is 
related to some problems studied by Darrell Duffie and coauthors (see Examples 2.1 and 3.3) 
in which they derive some models from the random local interactions between the financial 
agents involved, based on the ideas of DufRe and Sun (2007). Our initial goal was to provide 
a general framework in which this type of problems could be rigorously analyzed, and in 
particular prove a law of large numbers for them. In our general setting, W will be allowed 
to be any locally compact complete separable metric space. Considering type spaces of this 
generality is one of the main features of our model, and it allows us to provide a unified 
framework to deal with models of different nature (for instance, the model in Example 2.1 
has a finite type space and the limit solves a finite system of ordinary differential equations, 
while in Example 3.3 the type space is R and the limit solves a system of uncountably many 
integro-differential equations). 

To achieve this first goal, we based our model and techniques on ideas taken from the 
mathematical biology literature, and in particular on Fournier and Meleard (2004), where 
the authors study a model that describes a spatial ecological system where plants disperse 
seeds and die at rates that depend on the local population density, and obtain a deterministic 
limit similar to ours. We remark that, following their ideas, our results could be extended to 
systems with a non-constant population by adding assumptions which allow to control the 
growth of the population, but we have preferred to keep this part of the problem simple. 

The central limit result arose as a natural extension of this original question, but, as we 
already mentioned, it is much more delicate. The extra technical difficulties are related with 
the fact that the fluctuations of the process are signed measures (as opposed to the process 

which takes values in a space of probability measures), and the space of signed measures 
is not well suited for the study of convergence in distribution. The natural topology to 
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consider for this space in our setting, that of weak convergence, is in general not metrizable. 
One could try to regard this space as the Banach space dual of the space of continuous 
bounded functions on W and endow it with its operator norm, but this topology is too 
strong in general to obtain tightness for the fluctuations (observe that, in particular, the 
total mass of the fluctuations \/N(v^ — vt) is not a priori bounded uniformly in A^). To 
overcome this difficulty we will show convergence of the fluctuations as a process taking 
values in the dual of a suitable abstract Hilbcrt space of test functions. We will actually 
have to consider a sequence of cmbcddings of Banach and Hilbcrt spaces, which will help us 
in controlling the norm of the fluctuations. This approach is inspired by ideas introduced in 
Metivier (1987) to study weak convergence of some measure- valued processes using sequences 
of Sobolev cmbcddings. Our proof is based on Meleard (1998), where the author proves a 
similar central limit result for a system of interacting diffusions associated with Boltzmann 
equations. 

The rest of the paper is organized as follows. Section 2 contains the description of the 
general model. Section 3 presents the law of large numbers for our system, and Section 4 
presents the central limit theorem, together with the description of the extra assumptions 
and the functional analytical setting we will use to obtain it. All the proofs are contained in 
Section 5. 

2 Description of the Model 
2.1 Introductory example 

To introduce the basic features of our model and fix some ideas, we begin by presenting one 
of the basic examples we have in mind. 

Example 2.1. We consider the model for over-the-counter markets introduced in Duffie, 
Garleanu, and Pedersen (2005). There is a "consol", which is an asset paying dividends at 
a constant rate of 1, and there are N investors that can hold up to one unit of the asset. 
The total number of units of the asset remains constant in time, and the asset can be traded 
when the investors contact each other and when they are contacted by marketmakers. Each 
investor is characterized by whether he or she owns the asset or not, and by an intrinsic 
type that is "high" or "low" . Low-type investors have a holding cost when owning the asset, 
while high- type investors do not. These characteristics will be represented by the set of types 
W = {ho, hn, lo, In}, where h and I designate the high- and low-type of an investor while o 
and n designate whether an investor owns or not the asset. 

At some fixed rate A^, high- type investors change their type to low. This means that 
each investor runs a Poisson process with rate (independent from the others), and at 
each event of this process the investor changes his or her intrinsic type to low (nothing 
happens if the investor is already of low- type). Analogously, low-type investors change to 
high- type at some rate A^. The meetings between agents are defined as follows: each investor 
decides to look for another investor at rate (3 (understood as before, i.e., at the times of the 
events of a Poisson process with rate /?), chooses the investor uniformly among the set of A'' 
investors, and tries to trade. Additionally, each investor contacts a marketmaker at rate p. 
The marketmakers pair potential buyers and sellers, and the model assumes that this pairing 
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happens instantly. At equilibrium, the rate at which investors trade through marketmakers 
is p times the minimum between the fraction of investors willing to buy and the fraction of 
investors willing to sell (see Dufhe et al. (2005) for more details). In this model, the only 
encounters leading to a trade are those between hn- and Zo-agcnts, since high-type investors 
not owning the asset arc the only ones willing to buy, while low-type investors owning the 
asset are the only ones willing to sell. 

Theorem 1 will imply the following for this model: as N goes to infinity, the (random) 
evolution of the fraction of agents of each type converges to a deterministic limit which is 
the unique solution of the following system of ordinary differential equations: 



Here Uw{t) denotes the fraction of type-w investors at time t. This deterministic limit 
corresponds to the one proposed in Duffie et al. (2005) for this model (see the referred paper 
for the interpretation of this equations and more on this model). 

2.2 Description of the General Model 

We will denote by In = {l,---,-^} the set of particles in the system. In line with our 
original financial motivation, we will refer to these particles as the "agents" in the system 
(like the investors of the aforementioned example). The possible types for the agents will 
be represented by a locally compact Polish (i.e., separable, complete, metrizable) space W . 
Given a metric space ^{E) will denote the collection of probability measures on E, which 
will be endowed with the topology of weak convergence. When E = W, we will simply write 
V = V{W). We will denote by Va the subset of V consisting of purely atomic measures. 

The Markov process i^f we are interested in takes values in Va and describes the evolution 
of the distribution of the agents over the set of types. We recall that it is defined as 



where 5^, is the probability measure on W assigning mass 1 to w £ W and tj^ {i) corresponds 
to the type of the agent i at time t. In other words, the vector rj^ g W^'^ gives the 
configuration of the set of agents at time t, while for any Borel subset A of W, i'^ (A) is the 
fraction of agents whose type is in A at time t. 

The dynamics of the process is defined by the following rates: 

. Each agent decides to change its type at a certain rate j{w,v^) that depends on its 
current type w and the current distribution i/^ . The new type is chosen according to 
a probability measure a{w, , dw') on W . 

, Each agent contacts each other agent at a certain rate that depends on their current 
types wi and W2 and the current distribution i/^: the total rate at which a given type- 
wi agent contacts type-ii;2 agents is given by NX{wi,W2,i'^)i'^ {{w2})- After a pair 



(2.1) 



Uho{t) 
Uhn{t) 
Ulo{t) 
Uln{t) 



2(3uhn{t)uio{t) + pmm{uhn{t),Ulo{t)} + \uUlo[t) - XdUho[t), 
2j3uhn{t)uio{t) - pm:\n{uhn{t),Uio{t)} + XuUln{t) - XdUhnit), 
-2f3uhnit)uio{t) - pmm{uhn{t),Uio{t)} - XuUio{t) + XdUhoit), 
2j3uhn{t)uio{t) + pmin{uhn{t), Uio{t)} - XuUin{t) + XdUhn{t). 
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of agents meet, they choose together a new pair of types according to a probabihty 
measure b{wi,W2,i'^ ,dw'i^dw'2) (not necessarily symmetric in 'wi,W2) on WxW. 

For a fixed ^ G Va, a{w, ^, dw') and 6(^1, W2, dw['S^dw'2) can be interpreted, respectively, 

as the transition kernels of Markov chains in W and Wx W. 

Let B{W) be the collection of bounded measurable functions on W and Cb{W) be the 

collection of bounded continuous functions on W. For i/ G V and (p G B{W) (or, more 

generally, any measurable function ip) we write 



If dv. 



Observe that 



w 

N 



Wc make the following assumption: 
Assumption A. 

(Al) The rate functions ^{wjiy) and X{w,w',i') are defined for all v Cz V. They are non- 
negative, measurable in w and bounded respectively by constants 7 and A, and 
continuous in ly. 

(A2) a{w, ly, •) and b{w, w', z/, •) are measurable in w and w' . 

(A3) The mappings 

^{w,^) a{w,iy, ■) ^{dw) and 



w 



\{wi,W2,iy) b{wi,W2, ■) v{dw2) v{dwi), 



WJW 



which assign to each v £ Va & finite measure on W and W x W , respectively, are 
continuous with respect to the topology of weak convergence and Lipschitz with respect 
to the total variation norm: there are constants Ca,C'b > such that 



'y{w,iyi)a{'w,vi,-)iyi{d'w) - / ^{10,1^2)0(11], 1^2, ■)'^2{dw) 
w Jw 



< C'a||j^l-l'2||TV 



TV 



and 



X{wi,W2,iyi)b{wi,W2, Vl, ■) Vl{dW2) vi{dwi) 
WJW 



\{wi,W2, V2)b{wi,W2, V2,-) 1^2 (^^2) V2{dwi) 
WJW 



<Cb\\v^-iy2\\ 



TV- 



TV 



We recall that the total variation norm of a signed measure is defined by 

IIa'1Itv= sup 

f- llvlU<i 

(A3) is satisfied, in particular, whenever the rates do not depend on v. 
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3 Law of large numbers for 

Our first result shows that the process i'^ converges in distribution, as the number of agents 
N goes to infinity, to a deterministic limit that is characterized by a measure- valued system 
of differential equations (written in its weak form). 

Given a metric space S, we will denote by D{[Q,T], S) the space of cadlag functions 
V : [0, T] — > S, and we endow these spaces with the Skorohod topology (see Ethier and 
Kurtz (1986) or Billingsley (1999) for a reference on this topology and weak convergence 
in general). Observe that our processes i^^ have paths on D{[0,T],V) (recall that we are 
endowing V with the topology of weak convergence, which is metrizable). We will also denote 
by C([0,T], S) the space of continuous functions v: [0,T] — > S. 

Theorem 1. Suppose that Assumption A holds. For any given T > 0, consider the sequence 
of V -valued processes on [0,T], and assume that the sequence of initial distributions 
converges in distribution to some fixed I'o G V . Then the sequence converges in 
distribution in D{[Q,T]^V) to a deterministic vt in C([0,r],P), which is the unique solution 
of the following system of integro- differential equations: for every Lp g B{W) and t € [0,r], 

(i^t, V?) = (I'D, v?) + / / 'y{w,Us) / {(p{w') - Lp{w))a{w,iys,dw')i/s{dw)ds 
Jo Jw Jw 



(SI) 



/O JWJW JwxW 

■ b{wi, W2,Vs, dw[ ®dw'2) Vs{dw2) Vs{dwi) ds. 



Observe that, in particular, (SI) implies that for every Borel set A C W and almost every 

<G [o,r]. 



dvtiA) 



(^^{w,vt)+ J {\{w,w' ,vt) >^(w' iW,vt))vt{dw')^ vt(.dw) 



dt 

(SI') + / ^{w,vt)a{w,vt,A)vt{dw) 



w 

I X{w,w',i/t) b{w,w',i't,AxW) + b{w,w',i't,WxA) vt{dw') vt{dw). 
wJw 

Furthermore, standard measure theory arguments allow to show that the system (SI') ac- 
tually characterizes the solution of (SI) (by approximating the test functions (p in (SI) by 
simple functions). 

(SI') has an intuitive interpretation: the first term on the right side is the total rate at 
which agents leave the set of types A, the second term is the rate at which agents decide to 
change their types to a type in A, and the third term is the rate at which agents acquire 
types in A due to interactions between them. 

The following corollary of the previous result is useful when writing and analyzing the 
limiting equations (SI) or (SI') (sec, for instance. Example 3.3). 
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Corollary 3.1. In the context of Theorem 1, assume that vq is absolutely continuous with 
respect to some measure ji onW and that the measures 



I -f{w,i'o)a{w,VQ,-)voidw) and / / X{wi,W2,iyo)b{wi,W2,vo,-)vo{dwi)iyQ{dw2) 
Jw JwJw 

are absolutely continuous with respect to jj and pL® pL, respectively. Then the limit vt is 
absolutely continuous with respect to fi for all t £ [0,r]. 

The following two examples show two different kinds of models: one with a finite type 
space and the other with W = M.. The first model is the one given in Example 2.1. 

Example 3.2 (Continuation of Example 2.1). To translate into our framework the model 
for over-the-counter markets of Duffie et al. (2005), we take W = {ho, hn, lo, In} and consider 
a set of parameters 7, a, A, and b with all but A being independent of i'^ . Let 

^{ho) ^ "/(hn) ^ Xd, a{ho, ■) Sio, a{hn,-) = Sin, 

7(^0) = 7(/n) = A„, a{lo,-)^Sho, a{ln,-) = Shn- 

Observe that with this definition, high-type investors become low-type at rate A^ and low- 
type investors become high-type at rate A„, just as required. For the encounters between 
agents we take 

^ ^ [P if i^{{hn})iy{{lo}) ^ 0, 

b{hn,lo,i^,-) ^ 5/ho,in)^ and b{lo, hn, i^, ■) ^ dn„^hoi 



(where a A 6 = min{a, 6}), and for all other pairs wi,W2 S W, X{wi,W2, = (recall that 
the only encounters leading to a trade are those between hn- and Zo-agents and vice versa, 
in which case trade always occurs). The rates X{hn,lo,v) and X{lo,hn,v) have two terms: 
the rate (3 corresponding to the rate at which /in-agents contact /o-agents, plus a second 
rate reflecting trades carried out via a marketmaker. The form of this second rate assures 
that hn- and lo- agents meet through marketmakers at the right rate of pv{{hn}) A v{{lo}). 
It is not difficult to check that these parameters satisfy Assumption A, using the fact that 
X Ay ~ {x + y ~ \x ~ y\)/2 for x, y € K. 

Now let Uw{t) = vt{{w}), where vt is the limit of given by Theorem 1. We need to 
compute the right side of (SI') with A = {w} for each w e W . Take, for example, w ~ ho. 
We get 

Uho{t) = XuUlo{t)-XdUho{t)+liUhn{t)uio{t) + ^Uhn{t)f\Ulo{t)+[3uio{t)uh^^ 

which corresponds exactly to the first equation in (2.1). The other three equations follow 
similarly. 

Example 3.3. Our second example is based on the model for information percolation in 
large markets introduced in Duffie and Manso (2007). We will only describe the basic features 
of the model, for more details see the cited paper. There is a random variable X of concern to 
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all agents which has two possible values, "high" or "low" . Each agent holds some information 
about the outcome of X, and this information is summarized in a real number x which is 
a sufficient statistic for the posterior probability assigned by the agent (given his or her 
information) to the outcome of X being high. We take these statistics as the types of the 
agents (so W = M.). The model is set up so that these statistics satisfy the following: after 
a type-a;i agent and a type-X2 agent meet and share their information, xi + X2 becomes a 
sufficient statistic for the posterior distribution of X assigned by both agents given now their 
shared information. 

In this model the agents change types only after contacting other agents, so we take 
7 = 0, and encounters between agents occur at a constant rate A > 0. The transition kernel 
for the types of the agents after encounters is independent of and is given by 

b{xi,X2, ■) = b{x2,Xi,-) = (5(:ci+x2,xi+X2) 

for every Xi,X2 € M. This choice for the parameters trivially satisfies Assumption A. 

To compute the limit of the process, let A be a Borel subset of M. Then, since 7 = and 
A is constant, (SI') gives 



x-\-y ,x-\-y) 

(AxM)) vt{dy)vt{dx) 
= -2Xi^t{A) + 2\ [ S^+y{A) ut[dv) ut{dx) = -2\ut{A) + 2X [ vt{A - x) vt{dx), 



where A — x = \y ^9^: y + xG A\. Therefore, 

(3.1) vt{A) = -2\vt{A) + 2\{yt^vt){A). 

Using Corollary 3.1 we can write the last equation in a nicer form: if we assume that the 
initial condition vq is absolutely continuous with respect to the Lebesgue measure, then the 
measures vt have a density gt with respect to the Lebesgue measure, and we obtain 

/oc 
gt{z - x)gt{z)dz = -2\gt{x) + 2\{gt*gt){x). 
-00 

This is the system of integro-differential equations proposed in Duffie and Manso (2007) for 
this model (except for the factor of 2, which is omitted in that paper). 



4 Central limit theorem for z// 

Theorem 1 gives the law of large numbers for , in the sense that it obtains a deterministic 
limit for the process as the size of the market goes to infinity. We will see now that, under 
some additional hypotheses, we can also obtain a central limit result for our process: the 
fluctuations of v:^ around the limit vt are of order and they have, asymptotically, 

a Gaussian nature. As we mentioned in the Introduction, this result is much more delicate 
than Theorem 1. and we will need to work hard to find the right setting for it. 
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The fluctuations process is defined as follows: 

af = ^/]V(l.f -I.,). 

af is a sequence of finite signed measures, and our goal is to prove that it converges to the 
solution of a system of stochastic differential equations driven by a Gaussian process. As 
we explained in the Introduction, regarding the fluctuations process as taking values in the 
space of signed measures, and endowing this space with the topology of weak convergence 
(which corresponds to seeing the process as taking values in the Banach space dual of Cii{W) 
topologized with the weak* convergence) is not the right approach for this problem. The 
idea will be to replace the test function space Cb{W) by an appropriate Hilbert space Hi 
and regard cr^ as a linear functional acting on this space via the mapping ip i — > (c^, 93). In 
other words, we will regard as a process taking values in the dual TYi' of a Hilbert space 
Til. 

The space Hi that we choose will depend on the type space W. Actually, whenever W is 
not finite we will not need a single space, but a chain of seven spaces embedded in a certain 
structure. Our goal is to handle (at least) the following four possibilities for W: a finite set, 
Z'^, a "sufficiently smooth" compact subset of R'', and all of R'^. We wish to handle these cases 
under a unified framework, and this will require us to abstract the necessary assumptions on 
our seven spaces and the parameters of the model. We will do this in Sections 4.1 and 4.2, 
and then in Section 4.3 we will explain how to apply this abstract setting to the four type 
spaces W that we just mentioned. 

4.1 General setting 

During this and the next subsection we will assume as given the collection of spaces in which 
our problem will be embedded, and then wc will make some assumptions on the parameters 
of our process that will assure that they are compatible with the structure of these spaces. 
The idea of this part is that wc will try to impose as little as possible on these spaces, leaving 
their definition to be specified for the different cases of type space W. 
The elements we will use are the following: 

. Four separable Hilbert spaces of measurable functions on W, Hi, H2, H3, and 7^4. 

. Three Banach spaces of continuous functions on W, Co,C2, and C3. 

. Five continuous functions po, Pi, P2, Pa, P4 : W — > [l,oo) such that pi < pi+i for 
i = 0,1, 2, 3, Pi £ Ci for i = 0,2, 3, and for all w € W, p^iw) < Cp4{w) for some C > 
and p > 1 (this last requirement is very mild, as we will see in the examples below, but 
will be necessary in the proof of Theorem 2). 

The seven spaces and the five functions introduced above must be related in a specific 
way. First, wc assume that the following sequence of continuous embeddings holds: 

(Bl) Co ' — > Hi ' — > H2 ' — > C2 ' — > Hz = — > C3 ' — > 7^4, 

c 

where the c under the second arrow means that the embedding is compact. We recall that 
a continuous embedding Ei ^ E2 between two normed spaces Ei, E2 implies, in particular, 
that II • 11^:2 < C|| • 11^;^ for some C > 0, while saying that the embedding is compact means 
that every bounded set in Ei is compact in E2 ■ 
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Second, we assume that for i = 1, 2, 3, 4, if G Hi then 
(B2) |<^H| <C||^||„^p,H 

for all w e W, for some C > which docs not depend on ip. The same holds for the spaces 
d: for i = 0, 2, 3 and ip G Ci, 

(B3) \p{w)\<CMc^p,{w). 

The functions will typically appear as weighting functions in the definition of the norms 
of the spaces Tii and Ci. They will dictate the maximum growth rate allowed for functions 
in these spaces. 

We will denote by Hi' and C/ the topological duals of the spaces Hi and Ci, respec- 
tively, endowed with their operator norms (in particular, the spaces Hi and C/ are Hilbert 
and Banach spaces themselves). Observe that (Bl) implies the following dual continuous 
embeddings: 

(Bl') H4' ' — > C3' ' — > H3' ' — > C2 ' — > H2 ' — > Hi ' — > Co'. 

c 

Before continuing, let us describe briefly the main ideas behind the proof of our central 
limit theorem, which will help explain why this is a good setting for proving convergence of 
the fluctuations process. What we want to prove is that converges in distribution, as a 
process taking values in Hi , to the solution at of a certain stochastic differential equation 
(see (S2) below). The approach we will take to prove this (the proof of Theorem 1 follows 
an analogous line) is standard: we first prove that the sequence is tight, then we show 
that any limit point of this sequence satisfies the desired stochastic differential equation, 
and finally we prove that this equation has a unique solution (in distribution). To achieve 
this we will follow the line of proof of Meleard (1998). Our sequence of embeddings (Bl') 
corresponds there to a sequence of embeddings of weighted Sobolev spaces (see (3.11) in the 
cited paper); in particular, we will use a very similar sequence of spaces to deal with the case 
ly = K'^ in Section 4.3.4. One important difficulty with this approach is the following: the 
operator Js associated with the drift term of our limiting equation (see (4.1)), as well as the 
corresponding operators for (see (5.9)), cannot in general be taken to be bounded 
as operators acting on any of the spaces Hi- This forces us to introduce the spaces C,;, on 
which (B3) plus some assumptions on the rates of the process will assure that Js and are 
bounded. 

The scheme of proof will be roughly as follows. We will consider the semimartingale 
decomposition of the real- valued process (ct^,(^), for (f g Ha, and then show that the 
martingale part defines a martingale in 7^4 . This, together with a moment estimate on the 
norm of the martingale part in Hi and the boundedness of the operators in C3', will 
allow us to deduce that can be seen as a semimartingale in Tis', and moreover give its 
semimartingale decomposition. Next, we will give a uniform estimate (in N) of the norm of 

in C2 ■ This implies the same type of estimate in H2 , and this will allow us to obtain the 
tightness of in Hi . The fact that the embedding H2 ^ Hi is compact is crucial in this 
step. Then we will show that all limit points of cr^ have continuous paths in Hi and they all 
satisfy the desired stochastic differential equation (S2). Unfortunately, it will not be possible 
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to achieve this last part in Tii', due to the imboimdedness of Js in this space. Consequently, 
we are forced to embed the equation in the (bigger) space Co'. The boundedness of Js in Co' 
will also allow us to obtain uniqueness for the solutions of this equation in this space, thus 
finishing the proof. 

Our last assumption (D below) will assure that our abstract setting is compatible with 
the rates defining our process. Before that, we need to replace Assumptions (Al) and (A2) 
by stronger versions: 

Assumption C. 

(CI) There is a family of finite measures {T{w, z, ■)} ^ on W, whose total masses are 
bounded by 7, such that for every w £ W and every v £ V we have 

'y{w,i')a{w,i',dw') — / r{w, z,dw') i'{dz). 
Jw 

T{w, z, ■) is measurable in w and continuous in z. 

(C2) There is a family of measures {A(u'i, W2, z, ■)},^_^ ^^^r on WxW, whose total masses 
are bounded by A, such that for every wi,W2 S W and every v € V wc have 

X{'Wi,W2,i')b{wi,W2,v,dw[^dw'2) = / A{wi,W2, z,dw'i(>^dw'2) I'idz). 

J WxW 

A{wi, W2, z, •) is measurable in wi and 1^2 and continuous in z. 

The intuition behind this assumption is the following: the total rate at which a type-w 
agent becomes a type-w' agent is computed by averaging the effect that each agent in the 
market has on this rate for the given agent. Observe that, under this assumption, (A3) holds. 

Remark 4.1. Assumption C has the effect of linearizing the jump rates in i/. This turns 
out to be very convenient, because it will allow us to express the drift term of the stochastic 
differential equation describing the limiting fluctuations at ((S2) below) as JtCt for some 
Jt S Co' (see (4.1) and (5.9)). A more general approach would be to assume that the jump 
kernels, seen as operators acting on Co', are Frechet differentiable. In that case we would 
need to change the form of the drift operator Jt in the limiting equation and of Assumption 
D below, but the proof of Theorem 2 would still work, without any major modifications. To 
avoid extra complications, and since all the examples wc have in mind satisfy Assumption 
C, we will restrict ourselves to this simpler case. 

We introduce the following notation: given a measurable function tp on W, let 

rip{'w;z)= / {(p{w') — (p{w))T{w, z,dw') and 
Jw 

Aip{wi,W2; z) = / {ip{w[) + ip{w2) ~ tp{wi) - (p{w2)) A{wi,W2, z,dw[®dw'2). 

J WxW 

These quantities can be thought of as the jump kernels for the process associated with the 
effect of a type-z agent on the transition rates. Averaging these rates with respect to {dz) 
gives the total jump kernel for the process. 
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Assumption D. 

(Dl) There is a C > such that for aU w,z eW and i = 0,l, 2, 3, 4, 



p2 [uj') T{w,z,dw')<C {pj (w) + pj [z)) . 

w 

(D2) There is a C > such that for all wi,W2, z & W and i = 0, 1, 2, 3, 4, 

(p'K) + P?K)) K{wi,W2,z,dw[®dw'2) < C {pi{wi) + p^,{w2) + p1{z)) . 

WvW 

(D3) Let pi,p2 & V he such that < oo and define the foUowing operator acting on 

measurable functions (p on W: 

•^Mi,M2<P(-2^) = / Tfiw; z) p,i{dw) + / Tip{z]x) p2{dx) 

JW JW 

A(p{wi,W2; z) pi{dw2) pi{dwi) + / / Aip{w, z;x) pi{dw) p2idx) 

WJW JWJW 

A(p{z, w; x) p2 (dw) p2{dx). 

wJw 
Then: 

(i) Jfj.i,p,-> is a bounded operator on Ci, for i = 0,2,3. Moreover, its norm can be 
bounded imiformly in fii,p2- 

(ii) There is a C > such that given any ip ^ Co and any p-i, pi2, p-s, p4 € V satisfying 
{pi,pl) < oo, 

\\{Jt^uf^2 - •^M3,A'4)'/'llco - C'll'^llco (IIMi - A^sIIc^' + IIM2 - Ai4|lc2') • 

(Dl) and (D2) correspond to moment assumptions on the transition rates of the agents, 
and assure that the agents do not jump "too far". (D3.i) says two things: first, that the 
jump kernel defined by the rates preserves the structure of the spaces Ci and, second, that 
the resulting operator is bounded, which will imply the boundedness of the drift operators J, 
and mentioned above. (D3.ii) involves a sort of strengthening of the Lipschitz condition 
(A3) on the rates, and will be used to prove uniqueness for the limiting stochastic differential 
equation. Observe that by taking larger weighting functions pi, which corresponds to taking 
smaller spaces of test functions "Hi, we add more moment assumptions on the rates of the 
process; on the other hand, asking for more structure on the spaces Ti; and C^, such as 
differentiability in the Euclidean case, adds more requirements on the regularity of the rates. 

4.2 Statement of the theorem 

For ^ G "Hi (respectively Ci ) and <^ G H-i (respectively Ci) we will write 
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Given (p G Tii and z eW define 

Js(p{z)^ I Tip{w;z)vs{dw) + / Tip{z;x) Vs{dx) 
Jw Jw 



(4.1) + / / ^y^iwi,W2; z)iys{dw2)vs{d'Wi) 

JwJw 

[A^p{z, w; x) + Aiplw, z; x)] Vs{dw) Vs{dx) 

IWJW 

Observe that — Ju^Ms- Therefore, under moment assumptions on Vs, (D3.i) impUes that 
Js is a bounded operator on each of the spaces Ci. Observe that if we integrate the first 
and third terms on the right side of (4.1) with respect to Vs{dz), we obtain the integral term 
in (SI). In our central limit result, the variable z will be integrated against the limiting 
fluctuation process at. The other two terms in (4.1) correspond to fluctuations arising from 
the dependence of the rates on its other arguments (the types of the agents involved). 

The operator (or, more properly, its adjoint J*) will appear in the drift term of the 
stochastic differential equation describing the limiting fluctuations process, which will be 
expressed as a Bochner integral. We recall that these integrals are an extension of the 
Lebesgue integral to functions taking values on a Banach space, see Section V.5 in Yosida 
(1995) for details. 

Theorem 2. Assume that Assumptions C and D hold, that (Bl), (B2), and (B3) hold, and 
that 

, / , , , 2 \ 

< oo, 



y/N{v^ - J^o) ^ ^To, sup fif x/iV {v^ - ^^o) M 



sup E((i/^,p^)) < oo, and E((z^o,p^))< 



oo 



Ar>0 



hold, where the convergence in distribution above is in Jii . Then the sequence of processes 
converges in distribution in D{[0,T],'Hi') to a process at G C([0, T], Tii'). This process 
is the unique (in distribution) solution in Cq' of the following stochastic differential equation: 

(S2) at^ao+ / J*asds + Zt, 

Jo 

where the above is a Bochner integral, J* is the adjoint of the operator in Cq, and Zt is a 
centered Cq' -valued Gaussian process with quadratic covariations specified by 

[Z.{p,),Z.{ip2)\,= I I I [ {ipi{w')-^,{w)){ip2{w')-ip2{w))T{w,z,dw') 
Jo JwJwJw 

■ i's(dz) Vs{dw) ds 

+ ///// iipiiw[) + ipi{w2) ~ ifliwi) - ipi{w2)) 

Jo JwJwJwJwxw 

■ {f2{w[) + (^2(^2) - ¥'2(^1) - (^2(^2)) 

• A(wi,W2, z, dw'i (8) dw'2) Vs{dz) Vs{dw2) Vsidwi) ds 

for every ipi,ip2 ^ Co- 
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We will denote by Cf^'^'^ the sum of the two terms inside the time integrals above, so 

[Z.(^i),Z.(^2)],= fcr'^'^ds. 
Jo 

Remark 4.2. 

1. (S2) implies, in particular, that the solution at satisfies 

(S2-w) {crt,(p) = {ao,tp) + / (cts, Jsf) ds + Zt{tp) 

Jo 

simultaneously for every ip £ Cq. 

2. Observe that for any ipi, . . . ,ipk G Co, the process z!f^'"'''^^ = {Zt{ipi), . . . , Zt{ipk)) is 
a continuous M'^-valucd centered martingale with deterministic quadratic covariations, 
so it can be represented as 

Jo 

where [Cfj'^i'-''^* is the kxk matrix-valued process with entries given by \Cf'^'"'''^'']ij = 
^ ([C't]'^!' - '''''')^/^ is the square root of this matrix, and Bt is a standard k- 
dimensional Brownian motion. Thus, writing {at] i^i, . . . , tpk) — {{<Jt, 'fii) , ■ ■ ■ , {ct, 'fik)) 
we have 

(4.3) {at;cpi,...,^k) = f {at; Js^i, . . . , Js^fik) ds + ( {[C,]^'--'^>'f'^ dB,. 

Jo Jo 

3. The limiting fluctuations at have zero mass in the following sense: whenever 1 G Co 
and (cro,l) = 0, (ct*, 1) = for aU t G [0, T] almost surely. This follows from (4.3) 
simply by observing that, in this case, Jgl and are both always zero. 

Before presenting concrete examples where the setting and assumptions of this section 
hold, we present a general condition which allows to deduce that the assumptions (4.2) on 
the initial distributions ly^ , i/q, and aQ hold (namely, that Vq is a product measure). 

Theorem 3. In the setting of Theorem 2, assume that is the product of N copies of a 
fixed probability measure € V (i.e., is chosen by picking the initial type of each agent 
independently according to vq), and that¥.{{i/o, pf)) < oo. Then converges in distribution 
in V to vq, aQ converges in distribution in Hi to a centered Gaussian Tii -valued random 
variable a^, and all the assumptions in (4.2) are satisfied. 

4.3 Application to concrete type spaces 

In this part we will present conditions under which the assumptions of Theorem 2 are satisfied 
in the four cases discussed at the beginning of this section. 
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4.3.1 Finite W 



This is the easy case. The reason is that Ci,{W) can be identified with RI'^', and thus 
can be regarded as an Rl'^'l-vahied process, so most of the technical issues disappear. In 
particular, Theorem 2 can be proved in this case by arguments very similar to those leading 
to Theorem 1. 

In the abstract setting of Theorem 2, it is enough to choose pi = 1 and TLi = Ci ^ 
■j^\w\ ^ ("^(w^j for the right indices i £ {0, 1,2,3,4} in each case. (Bl) follows simply from 
the finite-dimensionality of RI'*'! and the equivalence of all norms in finite dimensions and 
(B2), (B3), and Assumption D are satisfied trivially. 

Theorem 2 takes a simpler form in this case. Write W — {wi, . . . , Wk}, 

k 

f^iit) = <^ti{wi\), Mcr) = Jsl{^o^}iwj) a-j, and 

i=i 

where a above is in M'"'. Also write F{(t) = (/i(<t), . . . , fk{<^)) and G{t) ~ {gij{t)) 
Observe that G{t) is a positive semidefinite matrix for all t > 0. 

Theorem 4a. In the above context, assume that Assumption C holds and that 




where the probability measures and Vt are taken here as elements of [0, 1]'^ and ao & MJ' . 
Then the sequence of processes cr^(t) converges in distribution in D{[0,T],M.'') to the unique 
solution cr(t) of the following system of stochastic differential equations: 

(S2-f) da{t) = F{<T{t)) dt + G^/^{t) dBt, 

where Bt is a standard k-dimensional Brownian motion. 

Example 4.3. This example provides a very simple model of agents changing their opinions 
on some issue of common interest, with rates of change depending on the "popularity" of each 
alternative. These opinions will be represented by = {— m, . . . , m} (m can be thought of 
as being the strongest agreement with some idea, as being neutral, and —m as being the 
strongest disagreement with it). Alternatively, one could think of the model as describing 
the locations of the agents, who move according to the density of agents at each site. 

The agents move in two ways. First, each agent feels attracted to other positions pro- 
portionally to the fraction of agents occupying them. Concretely, we assume that an agent 
at position i goes to position j at rate Pqijiy^ {{j}), where Q = {qij)ij^w is the transition 
matrix of a Markov chain on W. One interpretation of these rates is that each agent decides 
to try to change its position at rate /?, chooses a possible new position j according to Q, and 
then changes its position with probability vt{{j}) and stays put with probability 1 — J^t({j}). 
Second, each agent leaves its position at a rate proportional to the fraction of agents at its 
own position. We assume then that, in addition to the previous rates, each agent at posi- 
tion i goes to position j at rate apijiy^ {{i}), where P = {pij)ij^w is defined analogously 
to Q. This can be thought of as the agent leaving its position i due to "overcrowding" at 
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rate ai^tHi}) and choosing a new position according to P. We assume for simplicity that 
Pi,i = qi,i = for aU i € W. 

We win set up the rates using the notation of Assumption C: 



apij if k ^ i 

Pqij if k ^ j and A = 0. 

otherwise 



Assume that i/^ converges in distribution to some i^q E V, let vt be the limit given by 
Theorem 1 and write ut{i) ~ i^tdi}). It is easy to check that ut satisfies 



dufii) 



dt 

J — — m j — ~rn 



Now let at be the limit in distribution of the fluctuations process \/N\u^ — ut), and 
assume that the initial distributions i^q and vq satisfy the assumptions of Theorem 4a. It 
easy to check as before that 

m m 

Fi{(Tt) = 2a ^ pj.iUt{j)at{j) -2aut{i)at{i) + /3 ^ [qj^i - qij]{ut{i)at{j) + ut{j)at{i)). 

j—~7n j — — rn 

Thus, after computing the quadratic covariations we obtain the following: if ★ denotes the 
coordinate-wise product in RI^I (i-c, u * v{i) = u{i)v{i)) then the limiting fluctuations 
process at solves 

dat = 2aP^{ut*at) dt - 2aut*at dt + (3(^[Q^ - Q]crt)*ut dt 

+ f3(^[Q'-Q]ut)*atdt + ^/G(t) dBt , 

where Bt is a {2m + l)-dimensional standard Brownian motion and G{t) is given by 

{a'^Pk,iUt{k)'^ +aut{i)'^ + {qk,i + qi,k)ut{i)ut{k) if i = j 

~a{p.j^^ut{jf +p.i^jUt{if) ~ f3{qj^^ + q,j)ut{i)utij) if i ^ j. 



4.3.2 W ^Z'^ 

In this case Cb{W) is no longer finite-dimensional and, moreover, the type space is not 
compact, so we will need to make use of the weighting functions pi. We let D = [d/2\ + 1 
and take 

P,(.t) = ^l + \x\^'D. 
Clearly, we have in this case that < Cp4 for C — p — 2. 



16 



Consider the following spaces: 

Co = e°°{Z'^) = {ifi-.Z'^ such that \\ip\\^ < 00} , 
C, = i^'^'^iZ") = L: Z^^R such that = sup -^^^ < °° 1" = 2,3), 



2.iD 



= such that Hv'llLi? = ^ 1 + |^}L < ^} (* = 1,2,3,4) 



endowed with the norms defined within these definitions (we observe that pi does not appear 
explicitly in the definition of the spaces Ci, but the definition does not change if we replace 
the weighting function 1 + \xY^ appearing there by pi). These spaces are easily checked to 
be Banach (the d) and Hilbert (the Tii) as required. With these definitions we have the 
following continuous embcddings: 

(4.4) 

(these embeddings will be proved in the proof of Theorem 4b). 
To obtain (Dl) and (D2) wc will need to assume now that 

(4.5a) \ym^^ {y}) <C{1 + + l^D and 

(4.5b) Yl i\y^\''' + \y^-\''') A(xi,X2, z, {(2/i,2;2)}) < C {l + \x^r + 1^2^ + l^D 

for all xi,X2,z G Z'^ (the other six inequalities in (Dl) and (D2) follow from these two and 
Jensen's inequality). We remark that in Meleard (1998) the author also needs to assume 
moments of order 8D for the jump rates {8D + 2 in her case, sec (R'l) in her paper). 

Theorem 4b. In the above context, suppose that Assumption C holds and that (4.2), (4.5a), 
and (4.5b) hold. Then the conclusion of Theorem 2 is valid, i.e., converges in distribution 
in 7:>([0,r],£-2-D(Z'*)) (where ^-^■■^(Z'') is the dual of P^'^{Z'^)) to the unique solution at 
of the [Z'^y -valued) system given in (S2). 

We recall that the dual of £°°(Z'') can be identified with the space of finitely additive 
measures on Z'', and thus every <^ g £°°(Z'')' can be represented as {^i^)) .^.^j^d and we can 
write 

for if S £°°{Z'^). Therefore, (S2) can be expressed in this case in a manner analogous to 
(S2-f). 

Example 4.4. Here we consider a well-known model in mathematical biology, the Fleming- 
Viot process, which was originally introduced in Fleming and Viot (1979) as a stochastic 
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model in population genetics with a constant number of individuals which keeps track of the 
positions of the individuals. We will actually consider the version of this model studied in 
Ferrari and Marie (2007). 

Wc take as a type space W — Z+ and consider an infinite matrix Q = j))i.jgvFu{o} 
corresponding to the transition rates of a conservative continuous-time Markov process on 
VFU{0}, for which is an absorbing state (observe that, in particular, q{i, i) — — lih j))- 
Each individual moves independently according to Q, until it gets absorbed at 0. On absorp- 
tion, it chooses an individual uniformly from the population and jumps (instantaneously) 
to its position. We assume that the exit rates from each site are uniformly bounded, i.e., 
supj>i J2j£(wu{o})\{i} li^'j) < °° (t^^i^ ^° tli^^ (^1) satisfied). The rates take the fol- 
lowing form: 

{q{i,j) a j and i ^ j 

qihJ) + q(.i, 0) if fc = j and i^j and A = 0. 
ifi^j 

Observe that with this definition, the total rate at which a particle at i jumps to j when the 
whole population is at state v is given by q{i,j) + q{i, Q)v{{j}). 

We will write u^{i) = It is clear that this model satisfies the assumptions of 

Theorem 1. Therefore, if the initial distributions Uq converge, and wc denote by ut the limit 
given by Theorem 1, wc obtain that for each i > 1, 

^^^=^[g(z,j) + 9(*,0)u,(,)]^*(*). 

This limit was obtained in Theorem 1.2 of Ferrari and Marie (2007) (though there the 
convergence is proved for each fixed t). 

To study the fluctuations process we need to add the following moment assumption on 

Q: 

EA(*,j)<c(i+z«) 

for some C > independent of i. With this, if (4.2) holds, we can apply Theorem 2. By 
the remark following Theorem 4b, to express the limiting system for the fluctuations process 
it is enough to apply (S2-w) to functions of the form ip = li for each i > I. Doing this, 
and after some algebraic manipulations, we deduce that the limiting fluctuations process 
at is the unique process with paths in C([0, T], £°°(Z+)') satisfying the following stochastic 
differential equation: 

dat = Q'atdt+ {^Q{k,Q)(7t{k)\utdt + {^Q{k,Q)ut{k)\atdt+^/VtdBt, 

\k>l j \k>l j 

where Bt is an infinite vector of independent standard Brownian motions and Vt is given by 
[q{k,i) + q{k,0)ut{i)]ut{k) - [q{i,i) - q{i,0)]ut{i) + q{i,0)ut{i)'^ if i = j, 

Vt(,i,j) = { M» 

-q{hj)ut{i) - q{j,i)ut{j) - [g(i,0) + q{j,0)\ut{i)ut{j) if i ^ j. 
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4.3.3 W ^ fl, SL compact, sufficiently smooth subset of M 

Unlike the last case, the type space W is now compact, so we can simply take pi = 1. 
Nevertheless, W is not a discrete set now, and this leads us to use Sobolev spaces for our 
sequence of continuous embeddings: 

c 

(with D = ld/2\ + 1 as before), where C'^(O) is the space of continuous fimctions on fl with 
k continuous derivatives, endowed with the norm 

IMcin) = E sup|9"(p(a;)| , 

and H'^{Q) is the Sobolev space (with respect to the L^(il) norm) of order k, i.e., the space 
of functions on f2 with k weak derivatives in L^(fi), endowed with the norm 




The above embeddings are either direct or are consequences of the usual Sobolev embedding 
theorems, see Theorem 4.12 of Adams (2003). For these to hold we need to be sufficiently 
smooth (a locally Lipschitz boundary is enough). The compact embedding H^^{Vl) ^ 
H^{Q) is a consequence of the Rellich-Kondrakov Theorem (see Theorem 6.3 of Adams 
(2003)). 

In this case (Dl) and (D2) hold trivially. (D3) is much more delicate, and we will just 
leave it stated as it is. (The assumptions (H3), (H3)', and (H3)" of Meleard (1998) give 
some particular conditions which, if translated to our setting, would assure that (D3) holds. 
These conditions are suitable in her setting but they unfortunately rule out some interesting 
examples for us). 

Theorem 4c. In the above context, assume that Assumption and C holds, and that (D3) 
and (4.2) hold. Then the conclusion of Theorem 2 is valid, i.e., converges in distribution 
in £)([0,T],i7~3£'(f7)) (where H^^^{n) is the dual of H^'^iVi)) to the unique solution at of 
the (C^^ {^y -valued) system given in (S2). 

4.3.4 W 

This case combines both of the difficulties encountered before: W is neither discrete nor 
compact. To get around these problems we need to use now weighted Sobolev spaces. The 
weighting functions pi are given by 




where D = [d/2\ + 1 and q E N (to be chosen). We consider now the spaces C-'''^ of 
continuous functions ip with continuous partial derivatives up to order j and such that 
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lim|j.|^oo \d°'v{^)\ /(I + l^^l*^) = for all \a\ < j, with the norms 

\dMx)\ 



\a\<J- 

(as in Section 4.3.2, the weigthing functions pi do not appear explicitly here, but the definition 
does not change if we replace the term 1 + \x\'^ by -^/T+prp*^) and the weighted Sobolev 
spaces Wq'*^ (with respect to the norm) defined as follows: we define the norms 

n 1,2 V- f \dMx)f , 

\a\<3 

and let W^''^ be the closure in under this norm of the space of functions of class C°° with 
compact support. 

The right sequence of embeddings is now the following: 

Q3D,q , ^ yY^D,D+q ^ ^ T^2D.2D+q ^ ^ (jD.2D+q , ^ y^D,3D+q ^ ^ Q0,3D+q , ^ yyOAD+q 

g e N can be chosen depending on the specific example being analyzed: q = works for many 
examples, but as we will see in the next example, choosing a positive q {q = I in that case) 
can help, for instance, by making all constant functions be in C'^^''. These embeddings are, 
as before, either straightforward or consequences of the usual Sobolev embedding theorems 
(adapted now to the weighted case; see Meleard (1998), where the author uses the same type 
of embeddings, and see Kufner (1980) for a general discussion of weighted Sobolev spaces). 

To obtain (Dl) and (D2) we need to add the following moment assumptions on the rates, 
analogous to those we used in Theorem 4b: for all x,xi,X2, z G M'^, 

(4.6a) / |y|«^+2«r(x,z,dy)<C(l + |x|«^+2' + |z|«^+2«) and 



(4.6b) JR'ixR'' 

< C (1 + + |a;2|«^+2« + |z|«^+2«) . 

We observe that the power 8D + 2q appearing in this assumption corresponds exactly, when 
g = 1, to the moments of order 8Z? + 2 assumed in {H'l) in Meleard (1998). (D3), as in the 
previous case, is much more involved, so we will again leave it stated as it is. 

Theorem 4d. In the above context, assume moreover that Assumption C holds, and that 
(4.2), (D3), (4.6a), and (4.6b) hold. Then the conclusion of Theorem 2 is valid, i.e., 
converges in distribution in D(\Q,T],Wq '^^'^+'?) (where W^^^'^^'^ is the dual o/Wq^'^^"^) 
to the unique solution at of the ({C^^ •'^Y -valued) system given in (S2). 

Example 4.5 (Continuation of Example 3.3). In the previous section we obtained the system 
(3.1) that characterizes the information percolation model of Duffie and Manso (2007) by 
using (SI'). If we use (SI) instead we obtain 
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for all If g /B(M), where (ys^ip){z) = \^^^{x + z) Vsidx). 

To obtain the fluctuations limit, we need to check the assumptions of Theorem 4d. As 
we mentioned, we will take q ~ \. Assumption C holds trivially because \(w\^wi^v) and 
&(u'i, i(j2, i^, ■) do not depend on v. We will assume that the initial distribution of the system 
satisfies (4.2). (4.6a) and (4.6b) arc straightforward to check in this case. 

We are left checking (D3). Let E C"^'^ and take Mii /^2, /^3, G V having moments of 
order 10. We have that 

/oo />oo 
/ {2tp{wi + W2) - fiwi) - ip{w2)) fJ.i{dw2) Hi{dwi) 
-ocJ —00 

/oo 
{2ip{w + z) - ip{w) - (p{z)) [/ii(dw) + /i2(dw)] . 
-00 

The first term on the right side is constant in z, so it is in C^'^ (this is why we needed g = 1 in 
this example) . For the second term, since |(p(a;)| < C ||(y5||(,3.i (l + |a;|) and (/.t^, 1 + | • |-^*') < 00, 
the integral is bounded, and hence the derivatives with respect to z can be taken inside the 
integral, whence we get that this term is also in C^'^. The same argument can be repeated for 
and C°'4. The fact that the norm of this operator in these spaces is bounded uniformly 
in Hi, ^2 follows from the same argument. This gives (D3.i). Using the same formula it is 
easy to show that 

IIMI - M3||(C3-1)' + 11^2 - M4||(C3,1)' , 

which is stronger than (D3.ii). 

We have checked all the assumptions of Theorem 4d, so wc deduce that the fiuctuations 
process (j^ converges in distribution in W^^'^ to the unique solution of (S2) (which is an 
equation in (C"^'^)'). Writing down the formula for Jg in this case yields 

{(7s, Js'P) = 4 A {(Ts,t^s*<p) - 2 A {(7s,(fi) 

for every ip £ C"^'^. For the quadratic covariations we get 

^ (^jy^^ Vs*{fW2)) - 6A {vs, (^1) {vs, ip2) + 2A {vg, ^1^)2) 
for every (pi, (/32 € C'^'^. Therefore the limiting fluctuations satisfy 

(crt,<^) (cro,(/9) + A / [4(crs,:/s*'/') -2(crs,</j)]rfs + Zt((/3), 

with Zt being a centered Gaussian process taking values in the dual of C^'^ with quadratic 
covariations given by [Z{(pi), Z{(p2)]t = Jq Cf'^''^'^ ds for each tfi,tf2 € C^'^ ■ 

5 Proofs of the Results 

Throughout this section, C, Ci, and C2 will denote constants whose values might change 
from line to line. 
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5.1 Preliminary computations and proof of Theorem 1 

Since i/^ is a jump process in V with bounded jump rates, its generator is given by 

^Nf{i^)—N / j{w,i^) / ANf{iy;w;w')a{w,i',dw')iy{dw) 
Jw Jw 

^^'^^ +N / X{wi,W2,i') / ANf{iy;wi,W2;w[w'2) 

JwJw JWxW 

■ b{wi, W2, V, dwi ^dw'2)) v{d'Wi) h'(dw2) 

for any bounded measurable function / on where Ai\ff{i'; w; w') = f{v + N^^{6ii]' — 5wS) — 
f{v) and An f {i'\wi,W2;w[,w'2) = !{v + N^^{5u,[ + - - Sw^)) - f{v)- 

Given (f G B{W) we get by using (5.1) and Proposition IV. 1.7 of Ethier and Kurtz (1986) 
for f{y) = {v, Lp) that 

(:.f,^) = (z.o^,^)+Mf''^+ f [ 7(z^,^f) / {^{w')-v{w))a{w,,.^,dw')i.^{dw)ds 

Jo Jw Jw 

(^■2) + /7 f K^UW2,I^^) f {^{W[) + ^{W',) ^ V{W,) ~ ^{W2)) 

/o JwJw Jwxw 

b{wi,W2, ,dw'i®dw2) {dw2) {dwi) ds, 



where M^'''^ is a martingale starting at 0. This formula is the key to the proof of Theorem 1 
because, ignoring the martingale term, this equation has the exact form we need for obtaining 
(SI), and thus the idea will be to show that M^^''^ vanishes in the limit as TV ^ 00. This 
follows from the fact that the quadratic variation of '"^ is of order 0{1/N). More precisely, 
we have the following formula: for any ipi^ip2 S B{W), the predictable quadratic covariation 
between the martingales Mf. '"^^ and ''^'^ is given by 



■ a{w, , dw') (dw) ds 

(5.3) 1 



N 



A(wi,W2,I^s ) / {(pi{w[) + Lpiiw'^) - fl{wi) - Lpi{w2)) 

JWJW Jwxw 

■ {'P2{w[) + </?2(w2) - f2{wi) - (^2(^2)) 

• b{wi,'W2, , dw[ <S) dw'2) {d'W2) [dwi) ds. 

The proof of this formula is almost the same as that of Proposition 3.4 of Fournier and 
Meleard (2004) so we will omit it (there the computation is done for ipi = (^92, but the 
generalization is straightforward, and can also be obtained by polarization). 

Proof of Theorem 1. The proof is relatively standard, and its basic idea is the following. 
First one proves that the sequence of processes l^u^ ,ip) is tight in D([0,T],M) for each 
ip € Cf,(VF), which in turn implies the tightness of in £'([0, T], 7-"). The tightness of 
these processes follows from standard techniques and uses (5.2) and (5.3). Next, one uses 
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a martingale argument and (5.3) to show that any hmit point of i/^ satisfies (SI). FinaUy, 
using GronwaU's Lemma one deduces that (SI) has a unique solution. We refer the reader 
to the proof of Theorem 5.3 of Fournier and Meleard (2004) for the details. □ 

Proof of Corollary 3.1. Denote by {Tl^)i>o the sequence of stopping times corresponding to 
the jumps of the process i/^ . Let A be any Borel subset of W with fJ,{A) — and let (p be 
any positive function in B{W) whose support is contained in A. By (5.2), for every t G [0, T] 
we have that 

= E((^.o,^))+ErM^-f, 



(5.4) 




^{w,!^^ ) / {iflw') — (p{w)) a{w,i^g ,dw')iyg {dw)ds\ 
w Jw j 

A(wi,W2, i/f) / [ip(w\) + '.p{w'^) - ^}(wi) ~ ip{w2)) 
IWJW JwvW 

■ b{wi,W2, , dw'i ® dw'2 ) ) (dw2 ) {dwi ) ds^ . 

The first term on the right side of (5.4) is because the support of (p is contained A and 
i^oiA) = 0. The second term is by Doob's Optional Sampling Theorem. For the third term 
observe that for s < , — vq, so 



I {ip{w') - ip{w))a{w,vf ,dw')iyf{dw)ds 
w Jw 



E 



< 7E / / / \ip[w') - ip{w)\a{w,VQ,dw')vn{dw)\ 
\jo JwJw J 

which is since ip is supported inside A and the measure /^a(w', vq,-) h'o{dw') is absolutely 
continuous with respect to /i. The fourth term is by analogous reasons. We deduce that 
the expectation on the left side of (5.4) is 0, and therefore, since (p is positive, {v^^^n, p) = Q 



with probability 1. In particular, v^^^n is absolutely continuous with respect to /i for all 
t e [0, T] with probabihty 1. 

Using the strong Markov property we obtain inductively that {i^^^^n,^} = almost 
surely for every i > and t e [0,T]. Since the jump rates of the process are bounded, there 
are finitely many jumps before T with probability 1, and we deduce that {i'^ , (p) = almost 
surely for all t E [0, T]. Now if i/f is the limit in distribution of the sequence i/^ given by 
Theorem 1 and V3 G Cb{W), E{{iy^,ip)) {vt.ip) as iV ^ 00, so {vt, (^j) = for aU t e [0, T] 
whenever ip is supported inside A, and the result follows. □ 

5.2 Proof of Theorem 2 

We will assume throughout this part that all the assumptions of Theorem 2 hold. For 
simplicity, we will also assume that F = (these terms are easier to handle and are in fact 
a particular case of the ones corresponding to A). 
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Before getting started we recall that, by Parseval's identity, given any A £ Hi' and a 
complete orthonormal basis {(f>k)k>o of Hi, 

fc>0 

We will use this fact several times below. 



5.2.1 Moment estimates for and i^t 

Recall that we have assumed that sup^^Q E((i',5^, < oo and £((1^07^4)) < We need 
to show that these moment assumptions propagate to and i^t- 

Proposition 5.1. The following properties hold: 

(5.5a) sup E sup (y^ , < 00, and 

JV>0 YtG[o,T] ) 

(5.5b) sup {vt,(?^ < 00. 

te[o,T] 

The proof of this result will rely on an explicit construction of the process in terms of 
Poisson point measures. This is similar to what is done in Section 2.2 of Fournier and Meleard 
(2004) (though we will need to use a more abstract approach because our type spaces are 
not necessarily Euclidean), so we will only sketch the main ideas. 

We fix iV > and consider the following random objects, defined on a sufficiently large 
probability space: a T'-valucd random variable ly'Q (corresponding to the initial distribution) 
and a Poisson point measure Q{ds, di, dj, du, d9) on [0, Tjx/Arx/jvX [0, l]x [0, 1] with intensity 
measure (X/N) ds di dj du d9. We also consider a Blackwell-Dubins representation g oiV{Wx 
W) with respect to a uniform random variable on [0, 1], i.e., a continuous function g: ViWx 
W) X [0, 1] — > W xW such that g{£,, •) has distribution ^ (with respect to the Lebesgue 
measure on [0, f]) for all ^ G V{Wx W) and g{-,u) is continuous for almost every u E [0, 1] 
(see Blackwell and Dubins (1983) for the existence of such a function). This gives us an 
abstract way to use a uniform random variable to pick the pairs of types to which agents go 
after interacting. Finally, we introduce the following notation: rf {vP) will denote the i-th 
type, with respect to some fixed total order of W , appearing in (we recall that, under 
the axiom of choice, any set can be well-ordered, and hence totally ordered; moreover, this 
ordering can be taken to be measurable because W , being a Polish space, is measurably 
isomorphic to [0, 1]). With this definition, choosing a type uniformly from is the same as 
choosing i uniformly from and considering the type given by rii^{v^). Our process can be 
represented then as follows: 



,N _ ,,N 



JIn-jIn-jo Jo 



1 .1^ 

N 



^g^ (b(r)i(u^_),-ni (iyf_ ),i'f_ ,-),m) 

' ■'-6i<A(r,»(yf_),r,j(i/f_),i/f_)/A Q{ds,di,dj, du,d9), 
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where and are the first and second components of g (see Definition 2.5 in Fournier and 
Meleard (2004) for more details on this construction). 

Proof of Proposition 5.1. Since {v^ ^ pf) = {^^o , pf) + J2o<s<t [(^i^ ~ '^s^-^pT)] i it is easy to 
deduce from the last equation that 

{v^,pI) = {.^,pI)+ f f f C CUpi{g'm{.tu^{vti^^t,-\n)) 

Jo JIn-jIn-jo Jo '- 
+ pl(^.2(6(7?'(:.f_),r;^(z.f_),e, •),«)) -Pl('7^(e))-Pl(^^(^f-))" 
■ ■'■e<A(T,'(zyf_),TH(iyf_),i/f_)/A Q{ds,di,dj, du,d9) 
Taking expectations and ignoring the (positive) terms being subtracted we obtain 



EQup^^(.f,p^)^ <E((.o^,p^)) + ^ / E( ^^A(ry'(e),'7^(e),-. 



rp / N N 

[pl(ei(6(,f(i.f_),V(^f-),^f-, •),«)) 

+ pl{g\b{^\ut),if{vt).^f-.-).u)) 



du ds 



K{^o^pI))+I M I I I I [pl{w[) + pl{w'^)]A{w,,W2,z,dw[<E>dw'^) 
yJwJwJwJwxw 

■ v^{dz) v^{dw2) v^{dwi) I ds 



< 



K{{,y^,pl))+C n [pI{w,)+pI{w2) + pI{z)] 



IWJWJW 



iy';{dz)iy';idw2)iy';idwi) ds 



<E{{u^,pI))+C E sup {iyf,pi)]ds, 



where we used (D2) in the second inequality. By hypothesis E((i/(5^, pi)) is bounded imi- 
formly in N, so by Gronwall's Lemma we deduce that 

(5.6) e( sup {,^^,pI)] <Cie'^^^, 

\te[o,T] / 

with Ci and C2 being independent of N, whence (5.5a) follows. 

To get (5.5b), write {pi A L){w) = pKw) A L, and observe that, since pi A L e Cb{W), 

Theorem 1 imphcs that hmftr^oo IE(^suptg[Q (z^/^, pi A L)^ = sup^gjo.T] {'^t,pl A i), so by 
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(5.6) , 

sup {vt,p\KL) < Cie^'^. 
te[o,T] 

Using the Monotone Convergence Theorem it is easy to check that sup,g[Q (vg, p\ f\L) ^ 
supsg[o,T] {^siPi) as i ^ oo, and thus (5.5b) follows. □ 

For most of this section wc will continue ignoring the type-process 77^, working instead 
with the empirical distribution process wc arc interested in. However, wc will need 
to consider 77^ directly in Step 2 of the proof of Theorem 2, and we will need to use a 
moment estimate similar to (5.5a) for this process. Observe that statement of the theo- 
rem (and that of Theorem 1) makes no assumption on the distribution of rj^ , but instead 
only deals with the initial empirical distribution . Therefore we are free to choose rj^ 
in any way compatible with . For convenience we can construct 77;^ in the following 
way: assuming Vq takes a specific value £ Va, choose t/q (1) uniformly from Vq and 
then inductively choose 77^ (7) uniformly from the remaining — i -I- 1 individuals, i.e., 
from [NVq — <5^N(i) — • • • — /(^ — i + 1). It is clear then that, with this choice, 

77^ is exchangeable and jj^iLi^^l^ {i) = as required. Moreover, given any i G 
^{pHvo (*))) = '^(('^0^' pI)) J thus the moment assumption that we made on can be 
rewritten as supjv>o sup^g^^ E(^p1{riQ {i))) < 00 for all i G In- The proof of (5.5a) can then 
be adapted (by modifying slightly the explicit construction we made of ly^ to deal with 77^^) 
to obtain 

(5.7) sup sup E sup pliv^ii))] < 00. 

N>Oi£lN \te[o,T] J 

(We remark that the proof of this estimate uses (5.5a) itself). 
5.2.2 Extension of {v^ , •) and {i^t, ■) to Hi 

The 7^-valued process can be seen as a linear functional on B{W) via the mapping 
ip I — > {yP ^ip^, and the same can be done for vt- However, since 7^4 consists of measurable 
but not necessarily bounded functions, the integrals {y^ ^p)^ and {vt,'P>) may diverge. Our 
first task will be to show that these integrals are finite and, moreover, that (and vt) can 
be seen as taking values in 7^4' (and thus also in all the other dual spaces we are considering). 
A consequence of this will be that is well defined as an 7i4'-valued process. 

Proposition 5.2. The mapping if> e 7^4 ^ (yf is in H4' almost surely for every 
t € [0,T] and N > 0. Analogously, the mapping ip G 7^4 t-^ {i^t,^) is in TL^' for every 

[o,r]. 

Furthermore, Vt satisfies (SI) for every p) G 7^4, while satisfies (5.2) for every pi G 7^4 
almost surely. In particular, given any p G 7^4, Mf is a martingale starting at such that 
the predictable quadratic covariations (^M^''^'- , Al'^'^^'j^ are the ones given by the formula in 
(5.3) for all pi,(p2 G 7^4. 

Proof. We are only going to prove the assertions for , the ones for i/t can be checked 
similarly (and more easily). 
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The first claim follows directly from (B2) and Proposition 5.1: for ^p E Hi, 



Jw JW 

and the term inside the square root is almost surely bounded by (5.5a), so the mapping 
G 7^4 I — > iyt' 1 f) is continuous almost surely. 
Next we need to show that {y^ ,ip) satisfies (5.2) for all ip G TLi. That is. we need to 
show that the formula 



JWJWJW 



A(/3(wi, W2; z) {dz) {dw2) (dwi) ds 



N, 



.N I 



defines a martingale for each p G Ti^. Let (p G TL^ and to > and write ((/s A m){w) — 
tfiiw) A TO. tp /\ m is in B{W), so jVf/^''^^'" is a martingale. We deduce that given any 
< si < ■ ■ ■ < Sk < s < t and any continuous bounded functions ipi, . . . ,ipk on 7^4, if we let 



•••^.(-i! -Mf'^n 



then E(X'") = 0. Using the Monotone Convergence Theorem one can show that X' 

M^^J ■ ■ ■ M'^Ol^^t^'^ - ^i^'"^] as m ^ 00. On the other hand, the sequence (X™; 
is uniformly integrable. Indeed, using (B2) and (5.5a) one can show that 



m>0 



E 



(^4«)...V'fc«)[Mf''^^™-Mf''^-'"] )<Ct'Ei sup <oo. 



We deduce that 



E(Vi«) • • • Vfe«) [Mf ''^ - A^f '1 ) = lim mn = 0, 

-' rn— »oo 

which implies that M^'"^ is a martingale. The fact that ^M^''^% M^^'^^'j^ has the right form 
follows from the same arguments as those for (5.3) (here we need to replace (pi and (p2 by 
ip"^ and ip™ and then take to ^ 00 as above) . □ 



5.2.3 The drift term 

By Proposition 5.2, we have now that the fluctuations process is well defined as a process 
taking values in H4 and it satisfies 

+ Vn / / / / Atp{wi,W2;z) {dz)v^ {dw2)v^ [dwi) — Vs{dz)vs{dw2)i's{dwi) ds 
Jo JwJwJw L J 

for every ip G Ti^. The integral term can be rewritten as 
ft r r r 

Aip{wi,W2', z) (T^ {dz)v^ {dw2)i'^ {dwi) 



JwJwJw 



Vs{dz){(T^ {dw2)v^ {dwi) + Vs{dw2)(T^ (dwi)) 



ds. 
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Therefore, 

(5.8) (af,^) = V7V(.o"-z.o,^)+ViVMf--+ f\a^,J^^)ds 

Jo 

for each tp g 7^4, where 

(5.9) J^ip{z)^ ( ( Kip{wuW2;z)v^{d'W2)iy^{dwi) 

JwJw 

I Kip(w^ z]x) {dw) Vg{dx) + i I Atp{z,w,x) i'g(dw) ^^{dx). 
wJw ' JwJw 

Observe that = J^n ^i,^ and Jg = Jv,,^^ i where the operators J^i^^^^ ^'^^ the ones defined 
in Assumption D. Hence (D3) and Proposition 5.1 imply that and Js are bounded hnear 
operators on each space d {i = 0, 2, 3) and, moreover, for all ip G d, 

(5.10) lljf^ll^^ <C||^||c, and || J,,^!|c_ < C , 

almost surely for some constant C > independent of N and s. Similarly, given any if d Co, 

(5-11) II {J^ - Js) ^ ^ ll^llco \K - ^4c.' 

almost surely for some constant C > independent of N and s. 

5.2.4 Uniform estimate for the martingale term in Ti.4 

Proposition 5.2 implies that the martingale term M^'"*^ is well defined for all if e Ha. We 
will denote by M^^ the bounded linear functional on Ha given by [ip) = M^''^. 

Theorem 5.3. ^fNM^ is a cadlag square integrable martingale in Tii , whose Doob-Meyer 
process {{VNM^}^{ipi)){ip2) = N{y/NM^ {ipi),\/NM^ {ip2))^ (which is a linear operator 
from 7^4 to Ti.4 ) can be obtained from the formula in (5.3). Moreover, 



sup E 


sup 




2 ^ 


< OO 


Af>0 \ 


ytelo.T] 









Proof. We already know, by Proposition 5.2, that \/NAII^ is a martingale in H.^' with the 
right Doob-Meyer process. The fact that the paths of ^/NM^ are in D([0, T], 7^4') can be 
checked by the same arguments as those in the proof of Corollary 3.8 in Meleard (1998). So 
we only need to show the last assertion. Let {4'k)k>o be an orthonormal complete basis of 
7^4. We observe that, by (B2), if Xw G 7^4' is defined by Xwif) — </'(w) then 

J24iw) = \\x^,\\li,' <Cpliw). 

k>0 
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( sup 




\te[o,T] 





< E V sup N 
fc>o*e[o,T] 



= 4E 



JWJWJWJWyW 



Thus by Proposition 5.2 and Doob's inequality, 

fc>0 

A(wi, 1^2, z, dw'i ® dw'2) vf {dz) {dw2) [dwi) ds^ 

(P'K) + pI{w2) + P4K) + PIW2)) 

h-(wi, W2, z, dw'i ® dw'2) (dz) {dw2) [dwi)^ ds 
{2pl{wi) + 2pl{w2) + pI{z)) v^{dz) u^{dw2) i^^idwi)] ds 



< C 



E 



WJWJWJWxW 



< c 



E 



< C 



□ 



IWJWJW 

The last integral is bounded, uniformly in N, by Proposition 5.1. 
5.2.5 Evolution equation for in 'H3' 

Recall that our goal is to prove convergence of in D{[0,T],'Hi'). Therefore, a necessary 
previous step is to make sense of (5.8) as an equation in Hi - We will actually need to 
show something stronger: can be seen as a semimartingalc in Tis', whose semimartingale 
decomposition takes the form suggested by (5.8). We need the following simple result first 
(for its proof see Proposition 3.4 of Meleard (1998)): 

Lemma 5.4. For every N > there is a constant C{N) > such that 

sup E(||af||„,,) <C(iV). 
te[o,T] 

Recall that under our assumptions, need not be (and in general is not) a bounded 
operator on Ti.^, nor on any other Tii, and in fact (Tii) need not even be contained in 
Hi, so it does not make complete sense to speak of (J^) as the adjoint operator of . 
Nevertheless, for convenience we will abuse notation by writing {Jf)*cr^ to denote the 
linear functional defined by the following mapping: 



(Jf)'af(^) = (af,jf^)e 



Part of the proof of the following result will consist in showing that (jj^)* cr^ is actually in 



29 



Proposition 5.5. For each N > 0, is an H.^' -valued semimartingale, and its Dooh- 
Meyer decomposition is given by 

(5.12) = + + /Vf ds, 

Jo 

where the above is a Bochner integral in Ti^' . 

Proof. By Theorem 5.3 and the embedding Hi ^ , VnM^ is an 7^3 '-valued martingale. 
Thus, by (5.8), the only thing we need to show is that the integral term makes sense as a 
Bochner integral in 7^3'. The first step in doing this is to show that (Ji^) <J^ G Ti^' for all 
s G [0, T]. That is, we need to show that there is a C > such that 

(5.13) |(af,J»|<C||^||^^ 

for all ip e 7^3. Observe that by (5.10) and the embedding 7Y3 ^ C3, J^^ G C3 for e 7Y3, 
and thus 

|(af ,jf ^)| < ||af 11^^^, ||jf ^11^^^ <C\\a^\\^,^, 11^11^3 <C\\a^\\^,^, 

for such a function Lp by (Bl), so (5.13) holds almost surely by Lemma 5.4 and (Bl'). 

To see that the Bochner integral is (almost surely) well defined, we recall (see Section 
V.5 in Yosida (1995)) that it is enough to prove that: (i) given any function F in the dual 
of 7^3', the mapping s 1 — > F{i^J^Y a^) is measurable; and (ii) /g || (>^i^)*fi^||7^g/ ds < 00. 
(i) is satisfied by the continuity assumptions on the parameters and (ii) follows from (5.13), 
using the fact that the constant C there can be chosen uniformly in s. □ 

We omit the proof of the following corollary (sec Corollary 3.8 of Meleard (1998)): 

Corollary 5.6. For any N > 0, the process has paths in D{[0,T],'H3'). 

5.2.6 Uniform estimate for af on C2' 

Having given sense to equation (5.12) in 7^3', we can now give a uniform estimate for in 
C2' ■ This will be crucial for obtaining the tightness of in the proof of Theorem 2. 

Theorem 5.7. ^ 

sup sup E( llcrf'll ) < 00. 
N>OtG[0,T] ^ ' 

Proof. By (5.12) and the embedding 7^3' ^ C2 , 
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The first expectation on tlie right side is bounded uniformly in N by (4.2), and the same 
holds for the second one by (Bl') and Theorem 5.3. For the last expectation we have 























if 


/ tN\* n 
[Js ) 


ds 






Jo 




Jo 


C2' 





< T / E 



ds<CT I E sup ||cr 

\sG[0,t] 



^ \\C2 



dt. 



where we used Corollary V.5.1 of Yosida (1995) in the first inequality and (5.10) in the last 
one. Thus by Gronwall's Lemma we get E^supfgjQj^] llcrf^H^ ,^ < Cie^'^^ , uniformly in N, 
and the result follows. □ 



5.2.7 Proof of the theorem 

We are finally ready to prove Theorem 2. 

Proof of Theorem 2. As before, we will proceed in several steps. 

Step 1. Our first goal is to show that the sequence of processes is tight in D{[Q, T], Tii'). 
By Aldous' criterion (which we take from Theorem 2.2.2 in JofFe and Metivier (1986) and 
the corollary that precedes it in page 34) . we need to prove that the following two conditions 
hold: 

(tl) For every rational t G [0, T] and every e > 0, there is a compact K C Hi' such that 

sup F{af i K) < £. 

N>0 

(t2) If is the collection of stopping times with respect to the natural filtration associated 
to that arc almost surely bounded by T, then for every e > 0, 

lim limsup sup P( ||cr^^ — cr^ll^ , > e) = 0. 

r^O jv^oo s<r 

Observe that since the embedding of 7^2' into Jii is compact, (tl) will follow once we show 
that for any e > and t e [0,T] there is an L > such that sup^^qP^||(t^||^ , > < £■ 
This follows directly from Markov's inequality, (Bl'). and Theorem 5.7, since given any e > 0, 

sup P(||<Tf II >l)<^ sup E(||af II' ,) < -1 sup E(||af ||' ) < e 

N>0 ^ ^ ' ^ N>0 ^ 11^2/ L'' N>0 ^ 

for large enough L. 

To obtain (t2) we will use the semimartingale decomposition of in Tia' given in 
Proposition 5.5, i.e., = cr^ + VNM^^ + /J(jf )Vf ds. By Rebolledo's criterion (see 
Corollary 2.3.3 in Joffe and Metivier (1986)), (t2) is obtained for the martingale term ^/NM^ 
if it is proved for the trace of its Doob-Mcyer process {^/NM^}^ in Hi, and thus for erf if 

it is proved moreover for the finite variation term J^i^J^)* ds [a^ is tight by hypothesis). 
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We start with the martingale part. Let r be a stopping time bounded by T and let s > 0. 
Let (0fc)fc>o be an orthonormal complete basis of Using the same calculations as in the 
proof of Theorem 5.3 we get 



E 



< Cs, 



T + S 



X! {^k{w[) - <t)k{wi) + (t)k{w'2) - (t>k{w2)) 
WJWJWJWxW^yQ 

A(u;i, W2, z, dw'i ® dw'2) {dz) {d'W2) {dwi) 



uniformly in N . Thus by Markov's inequality. 



> 



< -Cs, 

e 



whence (t2) follows for the martingale term. 
For the integral term we have that 



Jo Jo 



dr 



< C 



■Hi' 
r+s 



E 



< E 



dr < Cs sup 

re[0,T] 



E Ik, 



N\ 



IIC2' 



for some C > 0, uniformly in N , where we used Corollary V.5.1 of Yosida (1995) as before 
and (Bl') in the first inequality and (5.10) in the second one. Using Markov's inequality as 
before and Theorem 5.7 we obtain (t2) for the integral term. 

Step 2. We have now that every subsequence of has a further subsequence which 
converges in distribution in -D([0, r],Hi'). Consider a convergent subsequence of erf, which 
we will still denote by erf, and let at be its limit in D([0, T], Hi'). Observe that the only 
jumps of fjf are those coming from and, with probability 1, at most two agents jump at 
the same time. Suppose that there is a jump at time involving agents i and j. Then given 

7^1, 

|(af ,^) - = |^(,7f (0) + V{V?{J)) V>iv^-{^)) V>iri^-{J))\ 



< 



C 



by (B2). We deduce by (5.7) that 
(5.15) e( 



\m\m 



II N 
sup (T" 

G[0,i] 



sup piiv^{i))+ sup pi(77f (.?■)) 

re[0,t] rG[0,t] 



N |K 
s- W-Hi' 



< 



c_ 

N 



and hence sup^^jQ \\a^ — cr^-W^ , converges in probability to as — > 00. Therefore, at 
is almost surely strongly continuous by Proposition 3.26 of Jacod and Shiryaev (1987). That 
is, we have shown that every limit point of erf is (almost surely) in C{[0,T],Ti.i). 
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Step 3. Our next goal is to prove that the sequence of martingales VNMjr converges in 
distribution in D{[Q,T],Ti.i') to the centered Gaussian process Zt defined in the statement 
of the theorem. That is, we need to show that given any ipi,ip2 G Ti-i, the sequence of 
R2-valued martingales //VAff '^'^'''^'^ = [yN Vn '"^^^ converges in distribution 

to {Zt{<p^),Zt{ip2))- 

By (5.12), ^/NM^^ and af have the same jumps, and thus (5.15) implies that 



(5.16) 



E sup 

\ se[o,t] 



0. 



On the other hand, wc claim that for every ipi,ip2 € Wi, 



(5.17) 



hm E 



ds. 



satisfies the hypotheses of the Martingale Cen- 



(5.16) and (5.17) imply that y/NMi 
tral Limit Theorem (sec Theorem VII. 1.4 in Ethier and Kurtz (1986)) so, assuming that (5.17) 

converges in distribution in -D([0, T],R^) to (Zt((^i), Zt{ip2))- 



holds, we get that ViVAff '^'^^''^^^ 

To prove (5.17) it is enough to consider the case ^pi ~ f2 ~ </5, the general case follows 
by polarization. Given ji G £'([0, T], Tii') let 



M^)= fill 

Jo JwJwJn 



WJWxW 



i(p{w[) + (^(^2) - <p{wi) - ip{w2)Y ^(wi, ii)2, dw\®dw'^) 

■ Hsidz) Hs{dw2) Hsidwi) ds. 



Then we need to prove that lim7v-»oo E(5't(i'^)) = '^t{v). Let p > 1 be the exponent we 
assumed to be such that < Cp4 for some C > 0.. Repeating the calculations in the proof 
of Theorem 5.3 and using Jensen's inequality we get that 

\^t{y'')\^ <\c^tMl,^ sup {vf,pl)Y<C2t^M\'^^ sup {u^^pI). 

Thus Proposition 5.1 implies that the sequence t{v^)) j^^^ is uniformly integrable, whence 
we deduce the desired convergence. 

Step 4. As in Step 2, let at be a limit point of . Observe that by the embedding 
Til ^ Co', (jf converges in distribution to at in D{[0, T],Co')- We want to prove now that 
at satisfies (S2-w). 

Fix ip e Co. By (5.12), 



{at,^) - {ao,Lp) - {as,Js'p)ds-Zt{(p) 



(5.18) 



TVMf - Ztiip)\ + [{at,^) - (af , ^)] + [{a^ , ^) - {a„,^)] 
+ A(af,J»-(af,J,^)]d.+ A(af,J,^)-(a„J,^)] ds 
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so we need to show that the right side converges in distribution to as iV ^ oo. The first 
term goes to by the previous step. The next two go to because at is a hmit point of 
and, since Js^ G Co, the last term goes to for the same reason. 

To show that the remaining term in (5.18) also goes to in distribution, it is enough to 
show that 



(5.19) 



E 



N- 



0. 



Since, by (5.10), — Js maps Co into itself, we get by using (Bl') and (5.11) that 



(.f , (Jf - .h) ^)\ < II (Jf - Js) < C ||af 11,^, ll^ll,^ 



.N 



^ II II II 



(5.19) now follows from this bound and Theorem 5.7. 

Step 5. We have shown in Step 4 that if at is any accumulation point of ct^, then at 
satisfies (S2-w) for every (y9 G Co. To see that the limit points of a^ actually solve (S2), the 
only thing left to show is that the integral term in (S2) makes sense as a Bochner integral in 
Co'. This can be verified by repeating the arguments of the proof of Proposition 5.5. 

Step 6. We want to prove now pathwise uniqueness for the solutions of (S2). Fix a centered 
Gaussian process Zt in Co' with the right covariancc structure and suppose that a] , a^ e Co 



are two solutions of (S2) for this choice of Zt 



Then a^ 



c^? = IoiJ>l ~ Jtf^l) ds, so 



sup 1 1 a 

te[o,T] 



ICo' 



< 



sup 

sG[0,t] 



I Co' 



dt. 



By (5.10), Js is a bounded operator on Co, and thus so is J* as an operator on Co'- Moreover, 
II J*||g^, can be bounded uniformly in s. Thus 



E sup \\al 
\te[o,T] 



ICo' 



< C 



E 



sup 

se[o,t] 



dt, 



and Gronwall's Lemma implies that al — erf for all t e [0,r] almost surely, so the pathwise 
uniqueness for (S2) follows. 

Step 7. We have now that any accumulation point at of the sequence af satisfies equation 
(S2), which has a unique pathwise solution. The last thing left to show is the uniqueness 
in law for the solutions of this equation. Since we have pathwise uniqueness, this can be 
obtained by adapting the Yamada-Watanabe Theorem to our setting (see Theorem IX. 1.7 of 
Revuz and Yor (1999)). The proof works in the same way assuming we can construct regular 
conditional probabilities in £)([0, T], Co'), which is possible in any complete metric space (see 
Theorem 1.4.12 of Durrett (1996)). This (together with the embedding Hi' ^ Co') implies 
that (S2) determines a unique process in C([0, T], Tii'). □ 
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5.3 Proof of Theorems 3 and 4a-4d 

Proof of Theorem 3. There are three conditions to check. The first one, =J> (Tq in Tii', 
follows directly from applying the Central Limit Theorem in R to each of the processes 
(a^ ,1^^ for e Til, while the condition sup^^g E(^i'(^, p^^) < oo is straightforward. For 

the remaining one we can prove something stronger, namely that sup^^Q E^||cr,^ ||^^,^ < oo. 

In fact, if (</'fe)fe>o is a complete orthonormal basis of 7^4 and 77,^ is chosen by picking the 
type r^Q (i) of each agent i G In independently according to i^o then 



E 



(j^O,0fe) 



A simple computation and (B2) (see the proof of Proposition 3.5 in Meleard (1998)) show 
that this is bounded by E((;^,^,p|)) + (vq^p^, which is in turn bounded by some C < 00 
uniformly in iV, so the result follows. □ 

For Theorems 4a (finite W), {W = Vl (ZW^ smooth and compact), and 4d {W = W^), 
we already explained why the assumptions of Theorem 2 hold, so the results follow directly 
from that theorem (together with (4.3) when W is finite). We are left with the case W = TJ^ . 

Proof of Theorem 4b. Let Lp e ^°°(Z'*). Then 

II ||2 \r- fix)'^ .^11 ||2 

Il'^ll2,15 = 2^ l+M^D < ^ll'/'lloo. 

where we used the fact that 2D > d implies that J^xei'^i^ + I^^P^)^^ < This gives 
the embedding £°°{'L'^) ^ (^-^{1/^). The other continuous embeddings in (4.4) are sim- 
ilar. To see that the embedding l'^-^{'L'^) ^ £2, 2D compact, observe that the family 
{ey)y^id C f^'^(Z'') defined by ey{x) = ^/l + \x\^lx=y defines an orthonormal complete 
basis of £'^'^{Z'^) and, using the same fact as above, 

SO the embedding is Hilbert-Schmidt, and hence compact. (B2) and (B3) follow directly 
from the definition of the spaces in this case. 

(Dl) and (D2) for pi are precisely what is assumed in Theorem 4b, and using this and 
Jensen's inequality we get the same estimates for pf, p2, and p^. We are left checking 
(D3). For simplicity we will assume here that A = 0. For (D3.i), the case Cq = £°°{Z'^) is 
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straightforward. Now if {fi,, 1 + | • \^^) < c», i = 1, 2, and if e ^°°-2^(Z'') 



1+ z 



\2D 



< 



+ 1 + U|2D E E (v'(2/)-V'W)r(z,x,{y})M2(W) 

cMU2d\^ E (i + l2/n r(x,z,M)(Mi({x}) + M2(W)) 



1 + |z 



2D 



< 



CM 



oo.2d 



1+ z 



2D 



1 + \zr + {1 + \xn 



,2D 



uniformly in z, where we used (4.5a) with a power of 2D instead of 8D. We deduce that 
||J"^i,M2<^lloo,2D ^ C'll'/^lloo.2D ^s required. The proof for i'^'^^(Z<^) is similar. For (D3.ii), 
consider £ i°°{Z'^) and /^i,M2,M3,M4 G V. Then 



w 



T ip{'w; z)\^lii{dw) — ^z{dw)\ + I Tip{z;w)\_^2{dw) — ^ii{dw)\ 



< - Ai3||foc,(2;d), + \\Tip{z; •)||^ \\^l2 - M4||^.=o(2 

Now ||r(y9(-; z)\\^ and ||r(y9(z; -jUj^ arc both bounded by 4A IjiysHoo, so we get 

as required. 



□ 
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